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ABSTRACT 


The  eqmvalence  of  the  finite  element  formulations  used  in  flexible  multibody  dynamics  is 
the  focus  of  this  investigation.  This  equivalence  will  be  used  to  address  several  fundamental 
issues  related  to  the  deformations,  flexible  body  coordinate  systems,  and  the  geometric 
centrifugal  stiffening  effect.  Two  conceptually  different  finite  element  formulations  that 
lead  to  exact  modeling  of  the  rigid  body  dynamics  will  be  used  in  this  investigation.  The 
first  one  is  the  absolute  nodal  coordinate  formulation,  in  which  beams  and  plates  can  be 
treated  as  isoparametric  elements.  This  formulation  leads  to  a  constant  and  symmetric  mass 
matrix  and  highly  nonlinear  elastic  forces.  It  is  demonstrated  in  this  study  that  different 
element  coordinate  systems  which  are  used  for  the  convenience  of  describing  the  element 
deformations  lead  to  similar  results  as  the  element  size  is  reduced.  In  particular,  two  element 
frames  are  used  in  this  study;  the  pinned  and  the  tangent  frames.  The  pinned  frame  has 
one  of  its  axes  passes  through  two  nodes  of  the  element,  while  the  tangent  frame  is  rigidly 
attached  to  one  of  the  ends  of  the  element.  Numerical  results  obtained  in  this  investigation 
using  these  two  different  frames  axe  found  to  be  in  a  good  agreement  as  the  element  size 
decreases.  The  relationship  between  the  coordinates  used  in  the  absolute  nodal  coordinate 
formulation  and  the  floating  frame  of  reference  formulation  is  presented.  This  relationship 
is  used  to  obtain  the  highly  nonlinear  expression  of  the  strain  energy  used  in  the  absolute 
nodal  coordinate  formulation  from  the  simple  energy  expression  used  in  the  floating  frame  of 
reference  formulation.  It  is  shown  in  this  paper  that  the  source  of  the  nonhnearity  is  due  to 
the  finite  rotation  of  the  element.  The  result  of  the  analysis  presented  in  this  paper  clearly 
demonstrates  that  the  instabihty  observed  in  high  speed  rotor  analytical  models  due  to  the 
neglect  of  the  geometric  centrifugal  stiffening  is  not  a  problem  inherent  to  a  particular  finite 
element  formulation.  Such  a  problem  can  only  be  avoided  by  considering  the  known  linear 
effect  of  the  geometric  centrifugal  stiffening  or  by  using  a  nonlinear  elastic  model  as  recently 
demonstrated.  Foiuier  analysis  of  the  solutions  obtained  in  this  investigation  also  sheds  new 
light  on  the  fundamental  problem  of  the  choice  of  the  deformable  body  coordinate  system 
in  the  floating  frame  of  reference  formulation. 
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1  INTRODUCTION 


The  floating  frame  of  reference  formulation  is  widely  used  in  flexible  multibody  simulations. 
In  this  formulation,  a  coordinate  system  is  assigned  to  each  flexible  body  in  the  multibody 
system  [13,  18].  The  location  and  orientation  of  this  body  coordinate  system  are  defined 
using  absolute  Cartesian  and  orientation  coordinates.  The  deformation  of  the  body  with 
respect  to  its  coordinate  system  is  defined  using  local  coordinates  that  can  be  introduced 
using  the  finite  element  method.  The  floating  frame  of  reference  formulation  leads  to  a  simple 
expression  for  the  elastic  forces  and  a  highly  nonlinear  expression  for  the  inertia  forces.  The 
simphcity  of  the  elastic  forces  is  due  to  the  fact  that  the  deformations  of  the  bodies  are 
defined  with  respect  to  their  respective  coordinate  systems.  The  nonlinearity  of  the  inertia 
forces  is  also  the  result  of  using  local  coordinates  to  describe  the  locations  of  the  points 
on  the  deformable  body.  In  order  to  define  the  absolute  position,  velocity  and  acceleration 
equations  in  the  floating  frame  of  reference  formulation,  a  coordinate  transformation  that 
defines  the  orientation  of  the  body  coordinate  systems  with  respect  to  a  selected  inertial 
frame  is  used.  The  use  of  this  coordinate  transformation  leads  to  the  nonlinearity  and  to  a 
strong  dynamic  coupling  between  the  absolute  reference  and  the  local  coordinates.  While  the 
finite  element  floating  frame  of  reference  formulation  can  lead  to  exact  modeling  of  the  rigid 
body  dynamics  when  beam  elements  are  used,  this  formulation  has  been  only  used  in  the 
analysis  of  small  deformation  problems  since  it  employes  non-isoparametric  beam  elements. 

Another  conceptually  different  formulation  which  is  introduced  recently  for  the  large  ro¬ 
tation  and  deformation  analysis  [2,  3, 10, 15, 19]  is  the  absolute  nodal  coordinate  formulation 
[8,  18].  In  this  formulation,  only  absolute  position  and  slope  coordinates  are  used.  Unlike 
the  floating  frame  of  reference  formulation,  the  absolute  nodal  coordinate  formulation  leads 
to  a  simple  expression  for  the  inertia  forces  and  a  highly  nonlinear  expression  for  the  elastic 
forces.  The  mass  matrix  is  constant  and  symmetric,  and  the  vector  of  centrifugal  and  Cori¬ 
olis  forces  is  identically  equal  to  zero.  In  the  absolute  nodal  coordinate  formulation,  exact 
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modeling  of  the  rigid  body  dynamics  can  be  obtained,  and  beams  and  plates  can  be  treated 
as  isoparametric  elements.  Furthermore,  the  formulation  of  joint  constraints  and  forces  can 
be  much  simpler  in  the  absolute  nodal  coordinate  formulation  as  compared  to  the  floating 
frame  of  reference  formulation  [8]. 

It  was  recently  demonstrated  that  the  highly  nonlinear  inertia  forces  obtained  in  the  float¬ 
ing  frame  of  reference  formulation  are  equivalent  to  the  simple  inertia  forces  obtained  using 
the  absolute  nodal  coordinate  formulation  [17].  The  result  of  this  study  clearly  demonstrated 
that  the  floating  frame  of  reference  formulation  does  not  lead  to  a  separation  between  the 
rigid  body  motion  and  the  elastic  deformation.  A  procedure  for  evaluating  the  inertia  shape 
integrals  required  to  evaluate  the  nonlinear  inertia  forces  of  the  floating  frame  of  reference 
formulation  from  the  constant  mass  matrix  of  the  absolute  nodal  coordinate  formulation  was 
presented  [17]. 

It  was  observed,  in  high  speed  rotor-craft  applications,  that  the  solutions  obtained  using 
the  floating  frame  of  reference  formiflation  incorrectly  exhibit  instabihty  problems  when  the 
angular  velocity  of  the  flexible  body  exceeds  a  certain  flmit  [11,  14].  This  incorrect  solution 
was  attributed  to  the  neglect  of  the  centrifugal  geometric  stiffening  effect.  Studies  which 
include  the  geometric  stiffening  effect  or  employ  a  nonlinear  elastic  model  showed  that  the 
instabihty  problem  can  be  solved  [7].  There  has  been  an  argument,  however,  that  the  in¬ 
stability  problem  is  an  inherent  problem  of  the  floating  frame  of  reference  formulation  and 
can  be  avoided  using  a  full  finite  element  representation;  and  in  such  a  full  finite  element 
representation,  a  correct  stable  solution  can  be  obtained  without  the  need  to  explicitly  ac¬ 
count  for  the  geometric  centrifugal  stiffening  effect  nor  the  elastic  nonlinearity.  The  analysis 
presented  in  this  paper  will  be  used  to  demonstrate  that  this  is  not  the  case.  To  this  end, 
the  absolute  nodal  coordinate  formulation  which  employs  a  full  finite  element  representation 
is  used.  Using  the  absolute  nodal  coordinate  formulation,  it  will  be  demonstrated  that  the 
generalized  elastic  forces  obtained  using  the  full  finite  element  representation  are  equiva- 
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lent  to  the  elastic  forces  obtained  using  the  floating  frame  of  reference  formulation,  thereby 
demonstrating  that  the  instability  problem  resulting  from  the  neglect  of  the  geometric  cen¬ 
trifugal  stiffening  effect  is  not  a  problem  inherent  to  only  the  floating  frame  of  reference 
formulation.  The  results  of  the  analysis  presented  in  this  paper,  with  the  results  previously 
reported  [17]  demonstrate  that  the  forces  used  in  the  floating  frame  of  reference  formulation 
can  be  obtained  from  the  forces  resulting  from  the  use  of  a  full  finite  element  representation 
by  simply  using  coordinate  transformation,  provided  that  no  linearization  of  the  kinematic 
equations  is  employed. 

This  paper  is  organized  as  follows.  In  Section  2,  a  brief  review  of  the  absolute  nodal 
coordinate  formulation  is  presented.  In  Section  3  and  4,  two  coordinate  systems,  the  pinned 
and  tangent  frames,  shown  in  Fig.  1,  are  used  to  define  the  element  deformations  and  the 
elastic  forces  in  the  absolute  nodal  coordinate  formulation.  In  Section  5,  it  is  shown  that 
the  results  obtained  using  these  two  frames  are  in  a  good  agreement  as  the  element  size 
decreases.  In  Section  6  and  7,  the  floating  frame  of  reference  formulation  is  briefly  discussed. 
In  Section  8,  the  relationship  between  the  coordinates  used  in  the  floating  frame  of  reference 
formulation  and  the  absolute  nodal  coordinate  formulation  is  presented.  This  relationship 
is  used  in  Section  9  to  demonstrate  the  equivalence  of  the  elastic  forces  used  in  the  floating 
frame  of  reference  formulation  and  the  absolute  nodal  coordinate  formulation.  In  Section  10, 
a  Fourier  analysis  of  the  solution  time  history  is  used  to  shed  more  light  on  the  fundamental 
problem  of  selecting  the  deformable  body  coordinate  system  in  the  floating  frame  of  reference 
formulation  [1,  4,  5,  9, 12, 16].  Summary  and  conclusions  drawn  from  this  study  are  presented 
in  Section  11. 
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2  ABSOLUTE  NODAL  COORDINATE  FORMULATION 


The  absolute  nodal  coordinate  formulation  is  a  non-incremental  finite  element  procedure  in 
which  the  coordinates  of  each  node  are  defined  in  a  fixed  inertial  coordinate  system;  hence  no 
transformation  matrices  are  required  to  define  the  kinematic  position  and  velocity  equations 
of  the  elements.  Using  the  motion  description  of  the  absolute  nodal  coordinate  formulation, 
a  constant  mass  matrix  is  obtained  for  the  finite  elements.  The  element  stiffness  matrix, 
however,  is  a  nonlinear  function  of  the  nodal  coordinates.  In  this  investigation  the  attention 
is  focused  on  beam  elements,  which  are  not  considered  as  isoparametric  in  the  classical  finite 
element  formulations  [10]. 

Referring  to  Fig.  2,  the  global  position  vector  r  of  an  arbitrary  point  P  on  the  element 
is  defined  in  terms  of  the  nodal  coordinates  and  the  element  shape  function  as 

=  Se.  (1) 

where  S  is  the  global  shape  function  which  has  a  complete  set  of  rigid  body  modes,  and  e  is 
the  vector  of  element  nodal  coordinates: 

T 

(2) 

This  vector  of  absolute  nodal  coordinates  includes  the  global  displacements 

Cl  =  ,  62  =  >  ^5  =  nL=i ,  ,  (3) 


Cl  62  63  64  65  66  67  63 


and  the  global  slopes  of  the  element  nodes,  that  are  defined  as 


dri 

dr2 

dri 

dr2 

J 

X— 0 

x=0 

x~l 

X=l 


(4) 


Here  x  is  the  coordinate  of  an  arbitrary  point  on  the  element  in  the  imdeformed  configuration, 
and  I  is  the  original  length  of  the  beam. 
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Since  absolute  coordinates  are  used,  a  cubic  polynomial  is  employed  to  describe  both 
components  of  the  displacements.  Therefore,  the  global  shape  function  S  can  be  written  as 

l-3e^  +  2f  0  +  0 

0  l-3f  +  2^3  0  + 

3^2  _2f  0  0 

0  3^2  _  2^3  Q 

where  ^  =  x/l.  It  can  be  shown  that  this  shape  function  contains  a  complete  set  of  rigid  body 
modes  that  can  describe  arbitrary  rigid  body  translational  and  rotational  displacements, 
provided  that  global  slope  coordinates  are  used  instead  of  infinitesimal  rotations.  Using  the 
absolute  coordinates  and  slopes,  it  can  also  be  shown  that  the  beam  element  defined  by  the 
shape  function  of  Eq.  (5)  is  an  isoparametric  element. 

Equations  of  motion  Using  the  global  position  vector  r  defined  in  Eq.  (1),  the  kinetic 
energy  of  the  finite  element  can  be  defined  as 

Iv  ~  (X  e  =  (6) 

where  p  and  V  are  the  mass  density  and  the  voliune  of  the  element  respectively,  and  Ma  is 
the  mass  matrix  defined  as  follows: 

Ma  =  /  (7) 

Jv 

This  is  a  constant  and  symmetric  matrix,  obtained  using  a  consistent  mass  approach.  It  can 
be  shown  [8]  that  this  matrix  leads  to  exact  modehng  of  the  rigid  body  inertia. 

In  order  to  develop  the  equations  of  motion  of  the  element,  the  vector  Qk  of  the  elastic 
forces  and  the  vector  Qa  of  the  externally  applied  forces  must  be  defined.  The  vector  Qfc 
can  be  defined  using  the  strain  energy  U  as  follows: 
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(8) 


As  it  will  be  demonstrated  later  in  this  section,  this  vector  is  a  highly  nonhnear  function  of 
the  vector  of  nodal  coordinates  e,  even  when  a  simple  linear  elastic  model  is  used. 

The  vector  Qa,  which  contains  the  generalized  external  forces,  including  the  gravity  force, 
can  be  defined  using  the  virtual  work  as 

8W,  =  Qj<5e.  (9) 

Using  the  expressions  of  the  kinetic  energy,  strain  energy,  and  the  virtual  work,  the  dynamic 
equations  of  the  finite  element  can  be  obtained  in  a  matrix  form  as  follows: 

Mae  +  Qk  =  Qa,  (10) 

where  all  the  terms  that  appear  in  this  equation  have  already  been  defined. 

Elastic  forces  In  order  to  demonstrate  the  nonlinearity  of  the  elastic  forces  even  when 
a  simple  hnear  elastic  model  is  used,  the  classical  Euler  Bernoulli  beam  theory  is  used  [20]. 
The  strain  energy  due  to  the  longitudinal  and  transverse  deformations  is  given  by 


where  E  is  Young’s  modulus,  A  is  the  cross  sectional  area,  /  is  the  second  moment  of  area, 
and  ui  and  Ut  are  respectively  the  longitudinal  and  transverse  deflections  as  shown  in  Fig. 
3. 

Let  i  and  j  be  unit  vectors  along  the  axes  of  the  beam  coordinate  system.  The  following 
relationship  holds  for  two-dimensional  problems: 

j  =  ii,  (12) 

where 

1  = 


0  -1 

1  0 


(13) 
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As  shown  in  Fig.  3,  the  location  of  an  arbitrary  point  P  on  the  beam  with  respect  to  the 
origin  O  of  the  local  frame  X1X2  is  defined  by 

u  =  rp  -  ro  =  (S  -  So)e,  (14) 


where  S  is  the  shape  matrix  given  by  Eq.  (5),  and  So  is  the  shape  matrix  evaluated  at  the 
reference  point  O.  The  longitudinal  and  transverse  defiections  of  point  P  can  be  defined  as 


Ul 

T* 

U^l  —  X 

u'^i  —  3;  1 

Ut  _ 

. . 1 

(15) 


According  to  these  definitions  of  the  longitudinal  and  transverse  deflections,  the  derivatives 
that  appear  in  the  strain  energy  expression  of  Eq.  (11)  are  defined  as  follows; 


and 


where 


T 

ii  =  (S"e)'^ii, 


S"  = 


d^s 

dx'^ 


By  defining  the  following  matrices 


Cl  =  t  S'^idx,  C2  =  t  S'^i  •  fS'dx,  C3  =  t  S"^Ii  •  i^V^S''dx, 
Jo  Jo  Jo 


(17) 


(18) 


(19) 


the  strain  energy  can  be  written  as 

U=^  [EA{1  -  2e'^Ci  +  e^Cse)  +  Ele^C^e]  .  (20) 

Clearly,  this  expression  is  not  a  quadratic  form  in  the  nodal  coordinates  vector  e,  as  in  the 
case  of  Unear  structural  mechanics.  Equation  (20)  is  a  highly  nonUnear  fxuiction  of  the  vector 
e  because  the  unit  vector  i  is  a  function  of  the  nodal  coordinates;  that  is  i  =  i(e).  The  form 
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of  the  matrices  given  in  Eq.  (19)  and  evaluated  using  the  shape  function  of  Eq.  (5)  are 
presented  in  the  Appendix. 

It  is  clear  from  the  analysis  presented  in  this  section  that  the  mass  matrix  that  results 
from  the  use  of  the  absolute  nodal  coordinate  formulation  takes  a  simple  form.  It  is  a 
constant  and  symmetric  matrix  ,  and  it  has  the  same  form  of  the  mass  matrix  used  in  hnear 
structural  dynamics.  The  elastic  forces,  on  the  other  hand,  have  a  much  more  complex 
expression  as  compared  to  the  vector  of  elastic  forces  used  in  linear  structmal  dynamics.  It 
is  one  of  the  objectives  of  this  investigation  to  establish  a  simple  procedure  for  evaluating 
the  nonlinear  elastic  forces  of  the  absolute  nodal  coordinate  formulation  from  the  simple 
expression  of  the  elastic  forces  used  in  linear  structural  dynamics.  To  this  end,  we  first 
demonstrate  that  different  element  coordinate  systems  used  to  formulate  the  strain  energy 
in  the  absolute  nodal  coordinate  formulation  lead  to  the  same  elastic  forces,  despite  the  fact 
that  the  definition  of  the  longitudinal  and  transverse  deformations  depend  on  the  choice  of  the 
coordinate  system.  As  a  consequence,  the  choice  of  the  element  coordinate  system  becomes 
immaterial  in  the  absolute  nodal  coordinate  formulation  as  the  element  size  decreases.  Using 
this  result,  we  demonstrate  in  later  sections  the  equivalence  of  the  elastic  forces  used  in  the 
absolute  nodal  coordinate  formulation  and  the  floating  frame  of  reference  formulation,  and 
use  this  equivalence  to  establish  a  simple  procedure  for  defining  the  nonlinear  elastic  forces 
of  the  absolute  nodal  coordinate  formulation  from  the  simple  expression  of  the  elastic  forces 
used  in  linear  structural  dynamics. 


3  THE  PINNED  FRAME 

In  order  to  describe  the  deformation  of  a  beam  element  in  the  absolute  nodal  coordinate 
formulation,  several  coordinate  systems  can  be  introduced.  It  is  important,  however,  to 
emphasize  that  the  use  of  a  coordinate  system  in  such  a  formulation  is  only  for  the  purpose 
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of  measuring  the  longitudinal  and  transverse  displacements  of  a  point  on  the  element  in 
order  to  evaluate  the  elastic  forces.  Therefore,  such  a  local  reference  frame  does  not  play 
the  same  role  as  in  the  case  of  the  floating  frame  of  reference  formulation.  For  instance, 
the  kinetic  energy  and  the  mass  matrix  used  in  the  absolute  nodal  coordinate  formulation 
were  formulated  in  the  preceding  section  without  the  need  for  any  local  coordinate  system. 
Furthermore,  it  is  theoretically  possible  to  develop  an  expression  for  the  strain  energy  and 
the  vector  of  elastic  forces  following  a  continuum  mechanics  approach  without  the  need  for 
introducing  a  local  reference  frame. 

In  the  floating  frame  of  reference  formulation  [18],  on  the  other  hand,  the  deformation 
of  the  bodies  is  often  described  using  elastic  modes,  and  the  elastic  coordinates  can  be 
considered  as  the  amplitudes  of  the  modes.  In  this  case,  where  low  frequency  modes  are 
often  used,  the  choice  of  the  boimdary  conditions  must  be  consistent  with  the  choice  of  the 
local  coordinate  system. 

In  this  section,  we  consider  the  case  of  a  pinned  frame,  which  has  one  of  its  axes  passes 
through  two  nodes  of  the  beam  element  as  shown  in  Fig.  lb.  The  longitudinal  and  transverse 
deformations  of  the  element  are  defined  in  this  frame,  and  they  are  used  to  define  the  elastic 
forces  in  the  absolute  nodal  coordinate  formulation.  In  the  following  section,  another  frame 
known  as  the  tangent  frame  is  employed.  The  tangent  frame  can  be  viewed  as  an  element- 
fixed  coordinate  system  which  has  one  of  its  axes  tangent  to  the  beam  at  one  of  the  nodes 
as  shown  in  Fig.  la.  In  Section  5,  a  comparison  between  the  pinned  and  the  tangent  frames 
is  presented,  and  it  is  shown  that  the  choice  of  such  frames  does  not  play  a  fundamental  role 
in  the  absolute  nodal  coordinate  formulation  as  in  the  case  of  the  floating  frame  of  reference 
formulation. 

The  pinned  frame  can  be  introduced  by  first  defining  the  unit  vector  i  along  a  line 
connecting  point  O  (the  origin  of  the  coordinate  system),  and  point  A.  Hence  the  imit 
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vector  i  is 


1  = 


ii 


r^-ro 


|rA-ro| 

Using  the  definition  of  the  vector  e  of  Eqs.  (2)-(4),  one  obtains 


(21) 


.  1 
^  d 


65 

66  —  62 

where  d  is  the  distance  between  the  nodes  at  0  and  A  defined  as 


(22) 


d=  y/{es  -  «i)^  +  (es  -  ^2)^- 


(23) 


The  derivatives  V’l  "02  of  H  and  i2  with  respect  to  the  nodal  coordinates  are 

■  (I)’ 


(24) 

(25) 


where  the  following  vectors  have  been  introduced 

lT 


Il  = 


-10001000 


,  h  = 


iT 


0-1000100 


(26) 


Substituting  Eqs.  (24)-(26)  into  Eq.  (8),  the  vector  of  elastic  forces  can  be  written  as 

Qfc  =  De-EA{hh+i2h)-EA{e'^li)'tP,-EA{e'^h)'il^2 

+{e^Die)iii>i  +  ^(e'’^D2e)(ii02  +  ^20i)  +  (e'^D3e)i202) 

where  the  matrices  Di,  D2,  D3  and  D  are  defined  as 


(27) 


Di=EACn+EIC3i,  =  EAC22  + EIC32,  D,  =  EAC23  +  EIC33,  (28) 


and 


D  =  (ii)2Di  +  «ii2D2  +  (»2)^D3. 


(29) 


The  matrices  are  defined  in  the  Appendix. 
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4  THE  TANGENT  FRAME 


In  this  section,  a  tangent  frame  that  can  be  viewed  as  an  element-fixed  coordinate  system  is 
employed.  This  frame,  shown  in  Fig.  3,  has  an  origin  rigidly  attached  to  the  node  point  0. 
In  this  case,  the  unit  vector  i  that  defines  the  Xi  axis  of  this  frame  is  given  by 


1  = 


es 

64 


(30) 


DiflFerentiation  of  the  direction  cosines  ii  and  i2  with  respect  to  the  vector  of  nodal  coordinates 
gives 

(31) 

(32) 


(dh 

V?1  = 

(  dh 

¥>2  = 

Ue 

T 


de  J  P 

T  , 


p 


-I4, 


-I4  -  -^l3, 

P 


where 


h  = 


0  0  1  0  0  0  0  0 


,  l4  = 


lT 


0  0  0  1  0  0  0  0 


(33) 


and 

/  =  +  {e,y  =  .  (34) 

The  quantity  y/r^r'  plays  the  role  of  the  longitudinal  deformation  gradient,  so  that  the  lon¬ 
gitudinal  strain  is  ei  =  —  1).  Hence  /  gives  a  measme  of  the  longitudinal  deformation 

that  occrus  at  the  left  end  of  the  beam,  at  point  O,  and  /  =  1  in  the  case  of  no  deformation. 
A  similar  comment  appfies  to  the  other  end  of  the  beam  element,  therefore  the  information 
provided  by  the  slopes  is  not  only  limited  to  information  on  the  direction  cosines. 

The  vector  of  elastic  forces  in  the  case  of  the  tangent  frame  is 


Qfc  =  E)e  —  d- 1212)  ~  EA{€p\i)ip^  —  EA[^^\2)'P2 

-t-(e'^Die)iiV?i  -t-  i(e'^D2e)(iiV?2  +  ^2^1)  +  (e^D3e)i2V52i  (35) 
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where  the  matrices  that  appear  in  this  equation  have  already  been  defined  by  Eqs.  (26), 
(28),  (29),  (31)  and  (32).  Note  that  this  expression  of  the  elastic  forces  becomes  identical  to 
Eq.  (27)  if  '01  and  02  replaced  by  and  respectively. 


5  COMPARISON  BETWEEN  THE  PINNED  FRAME  AND  THE  TANGENT 
FRAME 

Nmnerical  results  are  presented  in  this  section  in  order  to  compare  between  the  solutions 
obtained  using  the  pinned  and  tangent  frames.  The  study  model  used  is  the  four-bar  mech¬ 
anism  shown  in  Fig.  4.  The  mechanism  consists  of  a  crankshaft,  a  connecting  rod  (coupler), 
and  a  follower.  The  mechanism  is  driven  by  a  moment  applied  to  the  crankshaft.  The 
driving  moment  as  a  function  of  time  is  presented  in  Fig.  5.  The  inertia,  geometric,  and 
elastic  properties  of  the  components  of  the  four-bar  mechanism  are  shown  in  Table  1.  The 
table  shows  the  mass  m,  the  cross  sectional  area  A,  the  second  moment  of  area  I,  the  length 
I,  and  the  modulus  of  elasticity  E  of  the  mechanism  components. 

The  absolute  nodal  coordinate  formulation  is  used  to  obtain  the  nmnerical  results  in  the 
two  cases  of  the  pinned  and  tangent  frames.  Different  finite  element  models  that  employ 
different  numbers  of  elements  are  used  in  this  numerical  study.  Note  that  a  low  value  for  the 
modulus  of  elasticity  of  the  connecting  rod  is  used  in  order  to  allow  for  the  large  deformation, 
which  can  be  systematically  simulated  using  the  absolute  nodal  coordinate  formulation.  The 
results  of  the  computer  simulation  are  shown  in  Figs.  6-10.  In  aU  the  four  models  used, 
one  element  is  used  for  the  crankshaft  and  fom  elements  are  used  for  the  follower.  In  the 
first  model,  6  elements  are  used  for  the  coupler;  in  the  second  model,  twelve  elements  are 
used  for  the  coupler;  in  the  third  model,  eighteen  elements  are  used  for  the  coupler,  and  in 
the  fourth  model  thirty-six  elements  are  used  for  the  coupler.  Therefore,  the  first  model  has 
11  elements,  the  second  model  has  17  elements,  the  third  model  has  23  elements,  and  the 
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fourth  model  has  41  elements.  The  isoparametric  beam  elements  used  in  this  investigation 
are  defined  by  the  shape  function  of  Eq.  (5)  and  the  vector  of  nodal  coordinates  of  Eq.  (2). 
All  components  of  the  four-bar  mechanism  are  assumed  to  be  made  of  uniform  beams  which 
are  initially  straight.  The  gravity  effect  is  taken  into  consideration.  Pin  joints  are  used  to 
describe  the  connectivity  conditions  between  the  components  of  the  four-bar  mechanism.  In 
the  absolute  nodal  coordinate  formulation,  the  pin  joint  constraints  take  a  very  simple  form 
as  reported  by  Escalona  et  al  [8]. 

The  numerical  results  show  that  in  general  the  pinned  frame,  as  compared  to  the  tangent 
frame,  requires  a  less  number  of  elements  in  order  to  achieve  convergence.  For  instance, 
when  the  pinned  frame  is  used.  Model  3  in  which  the  coupler  is  divided  into  18  elements 
gives  results  which  are  in  a  good  agreement  with  the  results  obtained  by  discretizing  the 
coupler  into  36  finite  elements. 

Figure  6  shows  the  results  obtained  using  Model  1,  which  employs  the  least  ntunber  of 
elements  for  the  coupler.  These  results  show  that  there  is  a  significant  difference  between  the 
solutions  obtained  using  the  pinned  frame  and  the  tangent  frame  when  the  large  deformations 
of  the  coupler  are  considered. 

Figiures  7-9  show  the  effect  of  increasing  the  number  of  elements.  It  is  clear  from  the 
results  presented  in  these  three  figmes  that  as  the  size  of  the  elements  decreases,  better 
agreement  between  the  two  solutions  is  obtained.  In  particular,  the  results  shown  in  Fig.  9 
demonstrate  that  there  are  no  significant  differences  between  the  two  solutions  when  Model 
4  is  used. 

Figure  10  shows  the  transverse  deflection  of  the  midpoint  of  the  coupler  using  different 
models  and  different  frames.  The  transverse  deflection  is  determined  as  the  distance  of  the 
midpoint  to  the  straight  line  that  connects  the  ends  of  the  coupler.  It  is  interesting  to  note 
that  the  model  which  has  17  elements  (12  elements  for  the  coupler)  yields  good  results  when 
the  pinned  frame  is  used.  This  is  not  the  case  when  the  tangent  frame  is  considered.  More 
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elements  are  required  in  order  to  achieve  convergence.  The  good  convergence  characteristics 
of  the  pinned  frame  can  be  attributed  to  the  fact  that  the  deformation  within  the  element 
as  defined  with  respect  to  the  element  coordinate  system  remains  small,  despite  the  fact 
that  the  connecting  rod  rmdergoes  large  deformation.  This  is  not  the  case  when  the  tangent 
frame  is  used.  The  element  deformations  defined  with  respect  to  a  coordinate  system  fixed 
to  the  end  of  the  element  are  no  longer  small,  as  shown  by  Fig.  11.  A  computer  animation 
of  the  motion  of  the  four-bar  mechanism  is  presented  in  Fig.  12.  This  computer  animation 
is  obtained  using  Model  4  and  the  piimed  frame. 


6  FLOATING  FRAME  OF  REFERENCE  FORMULATION 

A  conceptually  different  formulation  from  the  absolute  nodal  coordinate  formulation  dis¬ 
cussed  in  the  preceding  sections  is  the  floating  frame  of  reference  formulation,  which  is 
widely  used  in  flexible  multibody  dynamics.  In  the  floating  frame  of  reference  formulation,  a 
mixed  set  of  absolute  and  local  deformation  coordinates  are  used  to  define  the  configuration 
of  the  deformable  body.  Each  body  in  the  multibody  system  is  assigned  a  body  coordinate 
system.  The  location  and  orientation  of  the  body  coordinate  system  in  an  inertial  frame  of 
reference  are  defined  using  a  set  of  reference  absolute  Cartesian  and  orientation  coordinates. 
The  deformation  of  the  body  with  respect  to  its  coordinate  system  can  be  described  using 
the  finite  element  method.  Since  the  element  shape  functions  include  rigid  body  modes,  and 
since  the  finite  displacements  of  the  body  are  described  using  the  reference  coordinates,  the 
rigid  body  modes  of  the  shape  function  must  be  eliminated  in  order  to  define  a  unique  dis¬ 
placement  field.  To  this  end,  a  set  of  reference  conditions  [18]  which  are  consistent  with  the 
kinematic  constraints  imposed  on  the  boundary  of  the  deformable  body  are  used.  Therefore, 
in  the  floating  frame  of  reference  formulation  the  choice  of  the  deformable  body  coordinate 
system  plays  a  fundamental  role  in  the  formulation  of  the  kinematic  position,  velocity  and 
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acceleration  equations.  Furthermore,  the  degree  of  complexity  of  the  inertia  forces  in  the 
floating  frame  of  reference  formulation  depends  on  the  choice  of  the  deformable  body  coor¬ 
dinate  system  [18].  This  is  not  the  case  when  the  absolute  nodal  coordinate  formulation  is 
used. 

In  general,  the  floating  frame  of  reference  formulation  leads  to  a  simple  expression  for 
the  elastic  forces  and  a  highly  nonlinear  expression  for  the  inertia  forces.  The  nonlinearity 
of  the  inertia  forces  is  due  to  the  coordinate  transformation  required  to  define  the  location 
of  an  arbitrary  point  on  the  body  in  the  global  inertial  frame  of  reference.  Such  a  coordinate 
transformation  is  not  required  in  the  absolute  nodal  coordinate  formulation,  and  as  a  con¬ 
sequence,  the  mass  matrix  in  the  absolute  nodal  coordinate  formrflation  takes  a  very  simple 
form.  Shabana  and  Schwertassek  [17]  demonstrated  the  equivalence  of  the  inertia  forces  used 
in  the  absolute  nodal  coordinate  formiflation  and  the  floating  frame  of  reference  formulation. 
It  was  shown  that  the  nonlinear  mass  matrix  of  the  floating  frame  of  reference  formulation 
can  be  systematically  obtained  from  the  constant  mass  matrix  of  the  absolute  nodal  coor¬ 
dinate  formulation.  Crucial  to  proving  this  equivalence  is  the  definition  of  the  coordinates 
in  the  floating  frame  of  reference  formiflation;  local  slopes  instead  of  infinitesimal  rotations 
must  be  used. 

The  fact  that  the  stiffness  matrix  in  the  floating  frame  of  reference  formulation  takes  a 
simple  form  leads  to  the  natural  and  interesting  question  of  whether  or  not  one  can  use  this 
simple  stiffness  matrix  to  formulate  the  nonlinear  elastic  forces  that  appear  in  the  absolute 
nodal  coordinate  formiflation.  In  order  to  answer  this  question,  one  needs  first  to  develop 
the  relationship  between  the  coordinates  used  in  the  absolute  nodal  coordinate  formulation 
and  the  floating  frame  of  reference  formulation.  Crucial  to  the  success  in  developing  this 
relationship  is  avoiding  the  kinematic  hnearization  resulting  from  the  use  of  the  infinitesimal 
rotations  as  local  deformation  coordinates  in  the  floating  frame  of  reference  formulation. 
Therefore,  in  the  following  section,  local  slopes  instead  of  infinitesimal  rotations  are  used 
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as  deformation  coordinates  to  formulate  the  strain  energy  in  the  floating  frame  of  reference 
formulation.  As  it  will  be  seen,  the  use  of  the  slopes  will  slightly  change  the  form  of  the 
strain  energy. 


7  STRAIN  ENERGY  IN  THE  FLOATING  FRAME  OF  REFERENCE  FOR¬ 
MULATION 


In  this  section,  an  expression  of  the  strain  energy  will  be  derived  for  beam  elements  using 
the  floating  frame  of  reference  formulation  in  terms  of  local  location  coordinates  instead 
of  the  familiar  expression  often  obtained  in  terms  of  deformation  coordinates.  This  new 
expression  will  be  used  in  the  following  sections  to  formulate  the  nonlinear  elastic  forces  in 
the  absolute  nodal  coordinate  formulation  from  the  simple  expression  of  the  elastic  forces 
used  in  the  floating  frame  of  reference  formulation.  To  this  end,  a  coordinate  system  that 
satisfies  the  simply  supported  end  conditions  is  considered.  The  location  and  orientation 
of  this  coordinate  system  is  identified  using  three  coordinates;  the  position  vector  R  of  the 
origin  O  and  the  angle  9  which  defines  the  orientation  of  the  beam  coordinate  system  with 
respect  to  a  global  inertial  system,  as  shown  in  Fig.  13.  These  three  coordinates  are  the 
reference  coordinates  which  can  be  written  as 


Qr  = 


R 

9 


It  follows  that  an  arbitrary  point  P,  whose  position  in  the  element  coordinate  system  is 
given  by  the  vector  u,  has  the  global  position  vector 


r  =  R  +  Au, 


(36) 


where  A  is  the  transformation  matrix 


A  = 


cos  9  —sin  9 
sin  9  cos  9 


(37) 
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In  the  case  of  a  rigid  body,  the  value  of  u  remains  constant  during  the  motion,  while  in 
the  case  of  deformable  bodies,  u  depends  on  time  and  can  be  expressed  in  terms  of  a  local 
shape-function  matrix  S/  and  the  vector  of  local  coordinates  q;  as  follows; 


u  = 


Ui 

U2 


=  S/qi. 


(38) 


In  order  to  be  consistent  with  the  assumed  displacement  field  used  previously  in  this  paper 
for  the  absolute  nodal  coordinate  formulation,  the  following  local  coordinates  are  employed 
for  the  beam  element  in  the  case  of  the  floating  frame  of  reference  formulation: 

r  iT 


qi 


91  92  93  94  95  > 


(39) 


where 


9i 


dui 

dx 


du2 

du2 

o’ 

X=0 

7 

a;=0 

-  /  ..N 

q3  =  Ui{A),  q4= 

X=l 

(40) 


in  which  point  A  is  the  right  node  of  the  beam  element  as  shown  in  Fig.  13.  Note  that  local 
slopes  instead  of  infinitesimal  rotations  are  employed,  in  order  to  avoid  the  linearization  of 
the  kinematic  equations.  The  local  shape  function  Si  is  given  by  (with  ^  =  x/l): 


0  m-2e+e)  o  o  i{e-e) 


(41) 


The  deflections  of  an  arbitrary  point  P  in  the  case  of  the  floating  frame  of  reference 
formulation  are  given  by 


" 

Ul 

= 

Ui  —  X 

Ut 

«2 

Using  this  equation,  the  expression  of  the  strain  energy  U f  (where  the  subscript  /  refers  to 
the  floating  frame  of  reference  formulation)  is  identical  to  Eq.  (11),  which  is  repeated  here 
for  convenience: 


(42) 
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In  this  case, 

(£)  -  wi 

(»)  -  |« 

where  Sn  and  S21  are  the  first  and  the  second  rows  of  the  local  shape  function  matrix  defined 
in  Eq.  (41),  and  the  prime  indicates  derivation  with  respect  to  x.  Substituting  the  preceding 
two  equations  into  Eq.  (42)  leads  to 

Uf  =  y‘  [EA  (1  -  2S;,q,  +  q?'sjs;,q,)  +  B/q^SS^S'^]  dx.  (45) 


Let 


T 


0  0  1  0  0, 


(46) 


2/2/15 

0 

-I/IO 

-/2/30 

0 

/•l  ,  EA 

0 

0 

0 

0 

0 

K,=EAlJ^  — 

-//lO 

0 

6/5 

-//lO 

0 

-/2/30 

0 

-//lO 

2/2/15 

0 

0 

0 

0 

0 

0 

(47) 


0  0  0  0  0 

0  4/^  0  0  2P 

0  0  0  0  0 

0  0  0  0  0 

0  2/^  0  0  4/2 


(48) 


Using  these  definitions,  the  expression  of  the  strain  energy  can  be  simply  written  as 


U,  =  ^EAl  -  Kjq,  +  iq,'^(Ki  +  K2)q,. 


(49) 


Note  that  the  strain  energy  contains  the  linear  term  ^EAl  —  K^qj.  Therefore,  the  strain 
energy  is  not  a  simple  quadratic  form  of  the  local  coordinates.  This  is  mainly  due  to  the  fact 
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that  the  local  coordinates,  as  defined  in  Eq.  (38),  describe  the  local  position  of  an  arbitrary- 
point  P  in  the  undeformed  state  as  well  as  the  deformation  of  this  point.  When  coordinates 
that  only  account  for  deformations  are  used,  the  simple  and  famiUar  quadratic  form  of  the 
strain  energy  can  be  obtained. 


8  RELATIONSHIP  BETWEEN  THE  COORDINATES 


In  this  section,  a  relationship  between  the  coordinates  used  in  the  absolute  nodal  coordinate 
formulation  and  the  floating  frame  of  reference  formulation  is  presented.  This  relationship 
will  be  used  in  the  following  section  to  obtain  the  nonlinear  elastic  forces  in  the  absolute 
nodal  coordinate  formulation  from  the  simple  expression  of  the  elastic  forces  in  the  floating 
frame  of  reference  approach. 

Using  Eq.  (40),  it  can  be  shown  that  the  nodal  coordinates  of  the  absolute  nodal  coor¬ 
dinate  formulation  can  be  expressed  in  terms  of  the  coordinates  used  in  the  floating  frame 
of  reference  formulation  as 


ei 

Ri 

62 

R2 

63 

qi  cos  9  —  q2  sin  9 

64 

qi  sin  9  +  q2  cos  9 

65 

Ri  +  93  cos  9 

66 

R2  +  93  sin  9 

67 

q4  cos  9  —  95  sin  9 

68 

94  sin  0  +  95  cos  9 

e(qr,q/)- 


(50) 
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Using  Fig.  14,  it  can  be  shown  that  the  inverse  relationship  is 

q-^  =  63  cos  9  +  sin  9 
q2  —  —63  sin  0  +  64  cos  9 

93  =  (es  -  ei)  cos 0  +  (eg  -  63)  sin 0  > ,  (51) 

94  =  67  cos  9  +  6$  sin  9 
q^  =  —67  sin  ^  +  eg  cos  9 

j 

where  the  angle  9  can  be  expressed  as  a  function  of  the  vector  e  using  Eq.  (22).  Equation 
(51)  can  be  written  in  a  matrix  form  as 

q/  =  Be,  (52) 

where 

0  0  cos  9  sm9  0  0  0  0 

0  0  —  sin0  cos  9  0  0  0  0 

B  =  —  cos  0  —  sin  ^  0  0  cos  9  sin  9  0  0  (53) 

0  0  0  0  0  0  cos  9  sin  9 

0  0  0  0  0  0— sin^  cos  9 

9  ELASTIC  FORCES 

It  is  shown  in  this  section,  using  the  coordinate  relationship  presented  in  the  preceding 
section,  that  the  nonlinear  elastic  forces  used  in  the  absolute  nodal  coordinate  formulation 
can  be  systematically  obtained  using  the  simple  expression  of  the  elastic  forces  used  in  the 
floating  frame  of  reference  formulation.  To  this  end,  Eq.  (52)  is  substituted  into  Eq.  (49) 
leading  to 

Uf  =  IeAI  -  KjBe  +  ^e'^B'^(Ki  +  K2)Be.  (54) 

Using  the  following  definitions: 

c  =  cos9,  s  =  sm9,  (55) 
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It  can  be  shown  that  this  expression  for  the  strain  energy  is  the  same  as  the  strain  energy  U 
of  Eq.  (20)  obtained  in  the  case  of  the  absolute  nodal  coordinate  formulation.  Considering 
the  expression  of  the  matrices  Ci,  C2,  and  C3  presented  in  the  Appendix,  it  can  be  shown 
that  Kio  =  EACi  and  Kj  =  EAC2.  However,  Kj  ^  EIC3.  On  the  other  hand,  the  simple 
structure  of  matrix  K*  is  consistent  with  the  type  of  reference  frame  used  for  the  evaluation 
of  the  longitudinal  and  transverse  deflections.  The  reference  frame  used  satisfies  the  simply- 
supported  end  conditions,  and  as  such,  a  variation  in  the  coordinates  of  the  end  points,  say 
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O  and  A,  makes  a  variation  in  the  rotation  9  of  the  frame,  but  does  not  affect  the  transverse 
deflections  which  depend  only  on  the  slopes.  This  explains  why  Kf  has  zeros  in  all  the 
positions  not  related  to  the  slopes. 

In  order  to  prove  that  the  expression  of  the  strain  energy  given  by  Eq.  (59)  is  the  same 
as  the  one  given  by  Eq.  (20),  it  is  enough  to  show  that  e^(Kt  —  EICs)e  is  identically  equal 
to  zero  for  every  e  and  for  every  length  1.  This  can  be  easily  demonstrated  by  carrying 
out  the  matrix  multipflcation,  and  noting  that  the  coordinates  of  node  A  can  be  written  as 
65  =  ei  4-  dcos^  and  6^  —  62  +  dsinO,  d  being  the  distance  between  the  nodes  as  shown  in 
Fig.  14. 

Straightforward  differentiation  of  Uf  with  respect  to  the  nodal  coordinates  gives  the 
expression  of  the  elastic  forces.  Note  that  instead  of  calculating  the  derivatives  of  K;o,  K; 
and  Kt  with  respect  to  6,  it  is  more  convenient  to  differentiate  with  respect  to  c  and  s,  where 
c  =  cos  9,  s  =  sin  9.  The  result  of  this  differentiation  is 

(^)  ^  “  (e'^K,o.c)V»i  -  (e'^K,o,.)V’2 

+Kze4-  ^(e'^K;,ce)t/>i  +  ^(e^K;,5e)V’2 

+Kfe  +  i(e'^Kt,ce)V>i  +  i(e^K(,5e)V’2,  (60) 

where 

i>.-»  ♦■■O' 

Here  V’l  and  tl^2  same  vectors  defined  by  Eqs.  (24)  and  (25),  since  c  =  ii  and  s  =  ^2. 

Note  that  Eq.  (60)  yields  the  same  results  as  Eq.  (27),  as  a  consequence  of  the  equivalence 
of  the  expressions  of  the  strain  energy  used  in  the  two  formulations. 

It  is  also  important  to  point  out  that  starting  with  Eq.  (50),  a  velocity  transformation 
matrix  can  be  defined  and  used  to  obtain  a  relation  similar  to  Eq.  (52),  but  with  a  matrix 
B  depending  on  q;  as  well  as  9.  This  approach  leads  to  a  more  complicated  expression  for 
the  elastic  forces.  On  the  other  hand,  the  approach  used  in  this  section  leads  to  a  matrix  B 
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that  depends  only  on  6. 


10  FREQUENCIES  AND  BOUNDARY  CONDITIONS 

The  equivalence  of  the  two  finite  element  formulations  discussed  in  the  preceding  sections 
can  be  used  to  shed  more  light  on  the  fimdamental  problem  of  selecting  the  deformable 
body  coordinate  system  in  the  floating  frame  of  reference  formulation.  In  the  absolute 
nodal  coordinate  formulation,  a  full  finite  element  representation  is  used.  In  this  case,  large 
deformation  problems  can  be  examined  as  previously  demonstrated  by  the  results  of  the 
four-bar  example  presented  in  this  paper.  The  floating  frame  of  reference  formulation,  on 
the  other  hand,  has  been  primarily  used  for  the  small  deformation  analysis.  This  is  due  to  the 
fact  that  non-isoparametric  beam  and  plate  elements,  which  employ  infinitesimal  rotations, 
have  been  often  used  with  the  floating  frame  of  reference  formulation.  Mode  reduction 
techniques  have  been  also  used  with  the  floating  frame  of  reference  formulation  in  order  to 
reduce  the  number  of  elastic  degrees  of  freedom. 

In  the  floating  frame  of  reference  formulation  interesting  and  fundamental  issues  related 
to  the  choice  of  the  boundary  conditions,  mode  shapes  and  deformable  body  coordinate 
systems  must  be  addressed  [16].  Different  sets  of  mode  shapes  that  correspond  to  different 
sets  of  end  conditions  and  natural  frequencies  can  be  chosen  for  a  flexible  link  to  yield 
approximately  the  same  results  as  previously  demonstrated  [16]. 

For  example,  the  deformation  of  the  connecting  rod  of  the  four-bar  mechanism  shown  in 
Fig.  4  can  be  modeled  using  mode  shapes  obtained  from  simply  supported  end  conditions, 
free-free  end  conditions,  or  double  cantilever  end  conditions.  The  first  six  mode  shapes  and 
the  corresponding  natural  frequencies  using  these  different  end  conditions  of  the  connecting 
rod  are  shown  in  Table  2  [6].  The  results  presented  in  Table  2  are  obtained  using  the 
dimensions  and  material  properties  of  the  coimecting  rod  reported  in  Section  5,  with  the 
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exception  of  the  moduli  of  elasticity.  For  the  connecting  rod,  E  is  increased  to  0.5E10 
MPa  in  order  to  obtain  small  deformations  and  justify  the  use  of  the  linear  model;  for  the 
crankshaft,  E  is  increased  to  I.OEIO  MPa;  and  for  the  follower,  E  is  increased  to  5.0E10 
MPa. 

As  demonstrated  in  previous  pubUcations,  the  simulation  results  obtained  using  these  end 
conditions  agree  well  [16].  Figmre  15  shows  a  comparison  between  the  results  obtained  using 
the  floating  frame  of  reference  formulation  and  the  absolute  nodal  coordinate  formrdation 
for  the  midpoint  deflection  of  the  coimecting  rod.  The  results  of  the  floating  frame  of 
reference  formulation  are  obtained  using  the  simply  supported  end  conditions  and  six  modes 
of  vibrations.  It  is  clear  from  this  figure  that  there  is  a  good  agreement  between  the  results 
obtained  using  the  two  formulations. 

With  this  good  agreement,  it  is  possible  to  examine  the  frequency  content  in  the  solution 
obtained  using  the  absolute  nodal  coordinate  formulation  which  does  not  employ  modes. 
These  frequency  contents  can  be  compared  with  the  natural  frequencies  presented  in  Table 
2  in  order  to  see  if  there  is  any  correlation  between  the  natural  frequencies  of  the  linear 
problem  and  the  frequency  contents  in  the  solution  of  the  nonlinear  multibody  problem. 
Figure  16  shows  the  Fourier  transform  of  the  solution  for  the  midpoint  deflection  of  the 
coimecting  rod.  Several  frequencies  appear  to  be  significant  from  the  results  presented  in 
this  figure.  It  is  difficult,  however,  to  draw  any  correlation  between  the  frequency  content 
in  the  solution  and  the  natural  frequencies  obtained  using  the  boundary  conditions  specified 
in  Table  2.  These  results  show  that  in  the  floating  frame  of  reference  formulation,  it  is 
important  to  choose  the  boundary  conditions  that  yield  a  shape  of  deformation  consistent 
with  the  kinematic  constraints  imposed  on  the  motion  of  the  deformable  body.  However, 
there  is  no  relationship  between  the  natural  frequencies  of  the  Unear  problem  that  result 
from  imposing  these  boimdary  conditions  and  the  frequency  content  in  the  solution  of  the 
nonlinear  problem. 
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11  SUMMARY  AND  CONCLUSIONS 


The  floating  frame  of  reference  formulation  is  widely  used  in  flexible  multibody  simulations. 
In  this  formulation,  a  coordinate  system  is  assigned  to  eax^h  flexible  body  in  the  multibody 
system.  The  location  and  orientation  of  this  body  coordinate  system  are  defined  using 
absolute  Cartesian  and  orientation  coordinates.  The  deformations  of  the  body  with  respect 
to  its  coordinate  system  are  defined  using  local  coordinates  that  can  be  introduced  using 
the  finite  element  method.  It  was  observed,  in  high  speed  rotor  craft  appfications,  that 
the  solutions  obtained  using  the  floating  frame  of  reference  formulation  incorrectly  exhibit 
instabihty  problems  when  the  angular  velocity  of  the  flexible  body  exceeds  a  certain  hmit. 
This  incorrect  solution  was  attributed  to  the  neglect  of  the  centrifugal  geometric  stiffening 
effect.  Studies  which  include  the  geometric  stiffening  effect  or  employ  a  nonlinear  elastic 
model  showed  that  the  instability  problem  can  be  solved.  There  has  been  an  argument, 
however,  that  the  instabihty  problem  is  a  characteristic  of  the  floating  frame  of  reference 
formulation  and  can  be  avoided  using  a  full  finite  element  representation  which  does  not 
account  for  the  geometric  centrifugal  stiffening  effect  nor  the  elastic  nonlinearity. 

The  analysis  presented  in  this  paper  demonstrates  otherwise.  To  this  end,  the  absolute 
nodal  coordinate  formulation  which  employs  a  full  finite  element  representation  is  used.  Un¬ 
like  the  floating  frame  of  reference  formulation,  in  the  absolute  nodal  coordinate  formulation 
global  position  and  slope  coordinates  axe  used.  The  formulation  leads  to  a  constant  mass 
matrix  and  a  highly  nonhnear  vector  of  elastic  forces.  Using  the  absolute  nodal  coordinate 
formulation,  it  was  demonstrated  that  the  generalized  elastic  forces  obtained  using  the  full 
finite  element  representation  are  equivalent  to  the  elastic  forces  obtained  using  the  fioating 
frame  of  reference  formulation,  thereby  demonstrating  that  the  instability  problem  resulting 
form  the  neglect  of  the  geometric  centrifugal  stiffening  effect  is  not  a  problem  inherent  to 
only  the  floating  frame  of  reference  formulation. 

The  results  of  the  analysis  presented  in  this  paper  with  the  results  previously  reported 
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[17]  clearly  demonstrate  that  all  the  forces  used  in  the  floating  frame  of  reference  formulation 
can  be  obtained  from  the  forces  resulting  from  the  use  of  a  full  finite  element  representation 
by  simply  using  coordinate  transformation,  provided  that  no  linearization  of  the  kinematic 
equations  is  employed.  Using  the  equivalence  between  the  absolute  nodal  coordinate  formu¬ 
lation  and  the  floating  frame  of  reference  formulation,  the  fundamental  problem  of  selecting 
the  deformable  body  coordinate  system  in  the  floating  frame  of  reference  formulation  is 
discussed. 
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APPENDIX 


The  stiffness  shape  integrals  that  appear  in  Eq.  (19)  are  as  follows; 


Cl  =  iili  +^2125  C2  =  (il)^C2l  +*1^2^22  +  (i2)^C23»  C3  —  (*l)^C 
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where 
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Table  1.  Parameters  used  in  the  simulation  of  the  four-bar  mechanism 


Body 

m  [kg] 

/[ml 

Z[m] 

E  [MPa] 

Crankshaft 

0.6811 

1.257E-03 

1.257E-07 

0.2 

l.OOOE+09 

Coupler 

2.4740 

1.960E-03 

3.068E-07 

0.9 

5.000E+06 

Follower 

1.4700 

7.068E-03 

3.976E-08 

0.5196174 

5.000E+08 

Table  2.  Mode  shapes  and  natural  frequencies  for  different  boundary  conditions 


Mode 

Simply  supported 

Shape  Freq.  [Hz] 

Double  cantilever 

Shape  Freq.  [Hz] 

Free-free^*^ 

Shape  Freq.  [Hz] 

1 

- \  45.81 

65.28 

^ - ^  103.85 

2 

183.24 

-  65.28 

286.26 

3 

412.30 

561.19 

4 

7^2  07 

927.67 

5 

1145.27 

^  ^  1145.50 

1385.78 

6 

/N  /N  IfidQ  10 

_  1145  50 

102SS1 

^  \y  \/  X/ 

The  rigid  body  mode  is  not  reported  in  this  table 


Fig.  6.  Large  deformations  of  the  four-bar  mechanism 
—  Model  1  with  pinned  frame 
—  Model  1  with  tangent  frame 
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Fig.  8.  Large  deformations  of  the  four-bar  mechanism 
—  Model  3  with  pinned  frame 
Model  3  with  tangent  frame 


Fig.  9.  Large  deformations  of  the  four-bar  mechanism 
—  Model  4  with  pinned  frame 
—  Model  4  with  tangent  frame 


Transverse  deflection  of  the  midpoint  of  the  connecting-rod  [m] 


a)  Pinned  frame 


b)  Tangent  frame 


Fig.  10.  Transverse  deformation  of  the  midpoint  of  the  connecting-rod  using  different  models 


Fig.  15.  Comparison  between  the  results  of  two  different  formulations 


